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ABSTRACT 



This thesis details the experimental analysis of an active 
damping control technique applied to the Naval Postgraduate 
School's Flexible Spacecraft Simulator using piezoceramic sen- 
sors and actuators. The mass property of the flexible arm is 
varied to study the frequency effects on the Positive Position 
Feedback (PPF) algorithm. Multi-modal dynamic response is 
analytically studied using a f inite- element model of a canti- 
levered beam while under the influence of three different 
control laws: a basic law derived from the Lyapunov Stability 
Theorem, PPF and Strain Rate Feedback (SRF) . The advantages 
and disadvantages of using PPF and SRF for active damping 
control are discussed. 
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I . INTRODUCTION 

A. BACKGROUND 

The evolution of spacecraft into larger and more complex 
structures, such as Space Station Freedom, has greatly expand- 
ed the problem of structural dynamics. The need to limit 
launch costs drives the engineer to a low mass design which 
results in these large structures being extremely flexible 
with very low frequency fundamental vibrational modes. These 
structural characteristics pose some difficult control prob- 
lems which are the focus of the emerging Control -Structures 
Interaction (CSIJ field. 

There are a variety of scenarios, such as slewing and 
pointing maneuvers, docking with other spacecraft and inter- 
actions with environmental forces that will produce unwanted 
structural vibrations. There may be requirements for sensor or 
antenna pointing or attitude control which can not be met with 
only the inherent structural damping in effect. This describes 
the need to devise techniques which actively control the dy- 
namic response of flexible structures. Any technique, in sim- 
plistic terms, is composed of two separate, but interrelated 
operations. First, the motion of the structure must be accur- 
ately measured by some type of sensor. Then, this information 
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must be fed back via an effective control algorithm to an ac- 
tuator on the structure that will counteract the vibrations. 

When trying to measure the dynamics of a structure, two 
possible quantities that can be measured are the acceleration 
of a point or points on the structural member and the strain 
within the member. However, with the low frequency vibrations 
characteristic of the structures we are interested in, precise 
acceleration measurements are difficult. This leaves strain as 
the quantity best suited for defining low frequency vibra- 
tions . 

There are a variety of sensors available for this applica- 
tion such as fiber optic sensors, strain gages and piezoelec- 
tric ceramics. Fiber optic sensors can be used to make 
measurements with a variety of methods, but are not yet prac- 
tical for complexity and technology problems. Strain gages are 
limited by a relatively low sensitivity and are highly subject 
to noise. The one type that has received much study and is 
seeing increased usage is the piezoelectric ceramic. Reference 
1 highlights its advantages as having: 

• high strain sensitivity 

• low noise 

• low to moderate temperature sensitivity 

• easy implementation 

The desired characteristics of a vibration suppression ac- 
tuator are low mass, infinite bandwidth, electrically powered 
and an internal force producer [Ref. 2] . A low mass 
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requirement is necessary so as not to significantly alter the 
physical characteristics of the structure as well as the 
overall launch mass. An ideal infinite bandwidth would theor- 
etically allow control of all vibrational modes of the struc- 
ture. Electric power reduces the complexity and mass inherent 
in a mechanical system. An actuator that produces internal 
forces only would greatly simplify the control algorithm for 
the entire structure. 

A number of materials have been developed for this purpose 
such as piezoelectric ceramic wafers (PZT) , electrorestrictive 
ceramic wafers (PMN) , piezoelectric polymer, film (PVF2) and 
shape memory metal (nitinol) wires. PZT and PMN can produce 
similar strain levels, but PMN is non-linear and very temper- 
ature dependent. Although nitinol can produce strains 300 
times greater than PZT, it is difficult to incorporate and not 
practical for most applications. When considering these trade- 
offs, PZT becomes the best choice for an actuator. Its suita- 
bility is summarized in the following characteristics [Ref .1] : 

• high stiffness 

• sufficient stress to control vibrations 

• good linearity 

• temperature insensitivity 

• easy to implement 

• low power consumption 
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To further investigate active damping control techniques, 
this technology has been incorporated into the Naval Postgrad- 
uate School's Flexible Spacecraft Simulator (FSS) . 

B. FOCUS OF THESIS 

The effectiveness and utility of using piezoelectric sen- 
sors and actuators to control vibrations will be analytically 
demonstrated via computer simulations. A model will be devel- 
oped for a simple cantilevered beam using the finite element 
method. The dynamic response of this system will be analyzed 
while under the influence of three different control laws: 

• a basic Lyapunov stability-derived control law 

• Positive Position Feedback (PPF) 

• Strain Rate Feedback (SRF) 

The Positive Position Feedback algorithm will also be in- 
vestigated experimentally. This is the particular control law 
which has been built into the Flexible Spacecraft Simulator. 
The mass property of the flexible arm can be adjusted to vary 
the natural frequency of the first vibrational mode. As will 
be shown, the effect of the piezoelectric actuator on the dy- 
namic response of the beam will change as the natural frequen- 
cy varies from a particular optimum value. This will highlight 
the advantages and disadvatages of this technique. 
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II. THEORETICAL ANALYSIS 



A. PIEZOELECTRIC THEORY 

Piezoelectric ceramics and crystals derive their utility 
in vibration control from the ability to convert electrical 
energy into mechanical energy and vice versa. When a force is 
applied to the material, the crystalline structure creates a 
voltage proportional to the force. Likewise, an electric po- 
tential applied across the material will alter its shape 
dimensionally. 

Man-made ceramics have an advantage over natural crystals 
in that they can be customized for each application. These 
ceramics can be manufactured in almost any size and shape and 
have their piezoelectric properties oriented in any specified 
direction. These properties are induced into the material dur- 
ing a process called "poling". While heated to a temperature 
above its Curie point, an electric field is applied across the 
material. This has the effect of realigning the crystalline 
structure such that it is elongated in a direction parallel to 
the electric field, the "poling" axis, and foreshortened in a 
direction perpendicular to the "poling" axis. Henceforth, as 
shown in Figure 1, any applied voltage of the same polarity as 
the "poling" voltage will elongate the ceramic along the 
"poling" axis and a voltage of opposite polarity will elongate 
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Figure 1 Piezoelectric actuator mode 



it in a direction perpendicular to the "poling" axis. In the 
reverse operation, a parallel compressive force or a perpen- 
dicular tensile force in relation to the "poling" axis will 
produce a voltage of the same polarity and vice versa. (See 
Figure 2 . ) [Ref . 3] 

The voltage (V s ) produced by a sensor is given by the 
relationship : 



- 



Ed, 3i 



t(e x + e 2 ) 



( 1 ) 



where E is the modulus of elasticity of the ceramic, d 3 , is the 
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Figure 2 Piezoelectric sensor mode 



lateral strain coefficient, D is the dielectric constant, t is 
the thickness and €, and e, are the lateral strains [Ref. 4]. 
This lays the groundwork for measuring the vibration of a 
structural member by converting the strain on the ceramic pro- 
duced by the flexural vibrations into a voltage which can be 
much more easily manipulated and processed. 

The effect of an actuator is described by the equation for 
the bending moment (m b (x,t)) at a cross-section of the struc- 
tural member: 
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m b (x, t) = Ely" (x, t) - k d v d ( t) 



( 2 ) 



where El is the flexural rigidity of the member, y"(x, t) is 
the second order partial derivative of the transverse dis- 
placement with respect to x (the longitudinal coordinate) , k d 
is a constant dependent on the ceramic piezoelectric constant, 
dielectric constant, width and the member thickness and v d (t) 
is the voltage applied to the ceramic [Ref. 5] . This shows how 
an applied voltage can produce a strain which, in effect, 
places concentrated moments (a couple) at the endpoints of the 
element on which the piezoceramic material is located. The 
direction of the, moment needed to counteract a vibration is 
changed by changing the polarity of the applied voltage. 

B. FINITE ELEMENT MODEL 

In order to effectively analyze the dynamic response of a 
system and the effect of an applied control law, an accurate 
and workable model of the system must be developed. The finite 
element method is used to discretize the equations of motion 
from unwieldy partial differential equations into second-order 
differential equations of generalized coordinates. 

1. Equations of Motion 

This analysis will consider a cantilevered beam as the 
system. From Bernoulli - Euler beam theory, the kinetic energy 
(T) , potential energy (U) , and work (W) done by an applied 
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moment from a piezoceramic actuator (M(x)) on each element 



are : 




(3) 




(4) 




(5) 



where x is the coordinate in the longitudinal direction, y is 
the coordinate in the transverse direction, V e is the trans- 
verse velocity, E is the modulus of elasticity, I is the mo- 
ment of inertia about the bending axis, p is the mass per unit 
length, Xj is the coordinate of the left end of the element and 
h is the length of the element. A typical finite element is 
shown in Figure 3 with the associated moment (M) and shear (V) 
forces . 

From the Extended Hamilton's Principle, the following 
relationship is used to formulate the equations of motion: 



Taking a variation of each term and assuming a constant moment 




( 6 ) 
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over the length of the element yields: 
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Figure 3 Finite Element Schematic 



Integrating the first term by parts gives the following varia- 
tional equation: 
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At this point, substitutions are made to discretize the motion 
into the form: 



y(x,t) -Ni(-x) q^t) (9) 

i=i 



The N f terms are third- order polynomials, chosen to satisfy the 
boundary conditions at both ends of the element, which define 
the mode shapes and are functions of x only [Ref. 6] : 




( 10 ) 



The q, terms are the system generalized coodinates which are 
functions of time only. The quantities they represent are: 

= Transverse Deflection at Left End of Element 
q 2 = Rotation at Left End of Element (11) 

g 3 = Transverse Deflection at Right End of Element 
g 4 = Rotation at Right End of Element 

Taking the appropriate derivatives of Equation (9) and substi- 
tuting back into Equation (8) and taking a summation of the 
ith terms yields: 
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f (ptf i (x)g i (fc)l^(x)5g.,(t) + EIN" (x) q d ( t)N'J(x) bq^ ( t ) ) dx 

J fc l Jx l 



- M J XJ+ n" (x)bq j (t) dx 

Jx l 



dfc = 0 , 



i=l 



( 12 ) 



Since this relationship must hold true over any time period 
and for any arbitrary variation 6q,-, the integrand must equal 
zero. Now, the familiar second- order equation of a vibratory 
system becomes evident in the form: 



My q ± (t) + K lj q i (t) = AGVj(x)| X f i , J = 1..J2 (13) 



The mass (M) and stiffness (K) matrices are computed from: 



My = p f^ b N 1 (x)N j (x) dx 
Ky = El [ Zl ' b n"(x)n"(x) dx 

J X, 



(14) 



Performing the integrations results in a system of four equa- 
tions for each element: 
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The complete mass and stiffness matrices for the entire system 
are obtained from a simple combination of the elemental matri- 
ces. Since the and q 4 coordinates of one element are the 
same as the q, and q 2 coordinates of the next element, the glo- 
bal matrix is formed by overlapping (adding) the upper left 
2x2 submatrix of an element with the lower right 2x2 submatrix 
of the previous element. Therefore, a n- element system has ma- 
trices with dimensions of (nx2)+2 by (nx2)+2. However, for a 
cantilevered beam, the q, and q 2 coordinates of the first ele- 
ment are always zero reducing the matrices to nx2 by nx2 . 

2 . Actuator Moment 



ment is derived using the Lyapunov Stability Theorem. The to- 
tal energy of the beam: 



is chosen as the Lyapunov function. This function meets the 
requirements of being positive definite, continuously 
differentiable and equal to zero at the equilibrium state. The 
only remaining condition needed to ensure stability is that 
the time derivative of U must be less than or equal to zero. 
[Ref. 7] 



The control moment produced by an actuator on an ele- 
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The time derivative of the energy function is: 
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( 17 ) 



Integrating this equation using the following relationships 
for a vibrating beam, bending moment (M) and shear (V) : 

p y + El = 0 (18) 

ax 4 



M = El 



dx 2 



(19) 



v . 4 ? - Ei^y 

dx dx 3 



( 20 ) 



while assuming no shear in the beam yields 



o-m 



( 21 ) 



The quantity in parentheses can be measured by a specially 
designed electric circuit which amplifies the current devel- 
oped over the sensor. To ensure that the derivative of the 
Lyapunov function is always less than zero, M must be of the 
form: 
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M = -k 



( 22 ) 



M is directly proportional to the applied voltage to the actu- 
ator such that k is a positive constant resulting in: 

& = -k £ 0 (23) 



In terms of the generalized coordinates, the moment is ex- 
pressed as: 

M = -k (g 4 - q 2 ) (24) 



The system of equations for an element with an actuator at- 
tached now becomes: 
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3 . Modal Coordinate Form 

The equations of motion can be decoupled and more eas- 
ily analyzed by transforming the generalized coordinates into 
modal coordinates. This transformation is of the form: 

q = $$ (26) 
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where is the weighted modal matrix and £ is a vector of 
modal coordinates. 

First, the eigenvalues and eigenvectors must be calcu- 
lated from the freely vibrating undamped system: 

Mg + Kg = 0 (27) 

For a n-dimensional system, there are n eigenvalues (X,-, i= 

l...n) and n eigenvectors (<£,., i=l...n) which satisfy the 
equation : 

- (28) 

The eigenvectors, also known as the normal modes, are orthog- 
onal with respect to the mass and stiffness matrices. This 
property allows the matrices to be diagonalized using the re- 
lationships : 

= M, (29) 

4>IJC<J) i = K l 



where the diagonal elements M,- and K ; , (i=l...n), are known as 
the generalized mass and stiffness elements for each mode. 

The weighted modal matrix (3>) is obtained by dividing 
the ith column of the modal matrix by the square root of the 
ith generalized mass: 
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Substituting the weighted modal matrix into Equation (29) pro- 
duces the diagonalizations : 



= I 
= K* 



(31) 



where I is the identity matrix and the diagonal elements of K* 
are the eigenvalues and squares of the modal frequencies. 

If the system is assumed to have viscous damping and 
excited by an arbitrary function Fu, the equations of motion 
take the form: 

Mq + Cq + Kq = Fu (32) 



If the damping matrix is proportional to either M or K, it is 
also diagonalized by the weighted modal matrix. Applying the 
transformation of Equation (26) to Equation (32) and premul- 
tiplying by $ results in n-uncoupled equations of the form: 
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C. COMPUTER SIMULATIONS 



To demonstrate the utility and stability of the specified 
control laws, the finite element model has been programmed 
using MATLAB . The basic approach is to put the modal equations 
into a state -space form and use MATLAB 's linear simulator 
function to produce a time response. The code for the PPF 
simulation is in the Appendix. The code for the other control 
laws is of a similar form except for the specific control law 
formulation . 

The following values for the physical parameters of the 
cantilevered beam were selected: 

• number of elements: n = 8 

• element length: h = .1 m 

• element density: p = .5 kg/m 3 

• product of modulus of elasticity and El = .9 Nm 2 

moment of inertia 

Using these values, the system mass and stiffness matrices 
were formulated and used to solve the eigenvalue problem. That 
solution provided the modal frequencies of the beam. The first 
three are: 

• C 0 i = 7.37 rad/s = 1.17 Hz 

• co 2 = 46.19 rad/s = 7.35 Hz 

• a> 3 = 129.42 rad/s = 20.6 Hz 

The inherent structural damping was introduced by arbitrarily 
damping ratios ( f ; ) such that the damping increased from the 
lowest value for the first mode to the highest value for the 
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.nth mode. In general, the high frequency modes are damped out 
more quickly then the lower frequency modes due to the inher- 
ent natural frequencies. 

For all the simulations, the beam was deflected and re- 
leased from rest such that the first three modes were excited. 
Simple geometric calculations were performed to determine the 
required initial values for each of the generalized coordi- 
nates. The open- loop response of the vibrating beam is shown 
in Figure 4 . 

From this point on, the simulations are customized to the 
specific control law utilized. The simulation development and 
results are detailed in the following sections for each of the 
three chosen control laws. 

1. Lyapunov Stability-Derived Control Law 

As shown in Chapter II. B. 2, the equations of motion 
for each element are of the form: 

+ K lj q 1 = Fk*q £ (34) 

where k* is a 1 x 4 vector containing the proportionality 
factor k: 

k*=[0k0-k] 05) 

F is a 4 x 1 vector dependent on actuator placement. If an 
actuator is located on an element, 
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II? MOTION - FEU OPEN LOOP 




TIME, sec 

Figure 4 Open- Loop Response 



F = [ 0 -1 0 1 ] T 



(36) 



Otherwise, it is a null vector. For the entire system of n 
elements, the individual 4x4 Fk* matrices are combined to 
form the n x n matrix F*. In modal form, the system is 
represented by: 

l + (D - F)t + = 0 (37) 



where : 
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The system is further manipulated into state space form: 

X = AX + Bu Y = CX + Du (3 8) 



where : 

0 J 
- K * F - D_ 



x - [ 5 1 r a = 



C is the identity matrix and B and D are null matrices. 

Four different simulations were run with varying actu- 
ator locations as follows: 

• Case 1 . Actuator at base (element l) 

• Case 2 . Actuator in middle (element 5) 

• Case 3 . Actuator at tip (element 8) 

• Case 4 . Combination of Cases 1-3 (elements 1,5 and 8) 

The results are shown in Figures 5-8. As expected, the actu- 
ator at the base has the most significant effect on damping 
vibrations since it can apply the largest moment. Although 
actuators in the middle and at the tip do provide significant 
damping when compared to the open- loop response, their contri- 
bution is minimal in the combined system. A comparison of Case 
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1 and Case 4 shows how much the base actuator dominates. This 



fact highlights the advantage of using a single actuator when 
considering cost and implementation issues. 



TIP DISPLACEMENT - TEM e=l 




TIME, sec 

Figure 5 Lyapunov / Actuator at base 



The proportionality factor k is a function of the 
physical constants of the piezoelectric which relate the ap- 
plied voltage to the moment created and a gain factor applied 
to the feedback circuit. For simplicity, it is assumed to be 
unity for these simulations. 

Varying k has a significant effect on the response. 
Theoretically, the system is stable for any positive value. 
However, larger values are, in effect, placing a moment with 
a magnitude larger than required which tends to "overcontrol" 
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TIP DISPLACEMENT - FEU k = l e=S 




Figure 6 Lyapunov / Actuator in middle 



TIP DISPLACEMENT - FEM krl e=8 




Figure 7 Lyapunov / Actuator at tip 
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TIP DISPLACEMENT - PEN k=l e=1.5.8 




Figure 8 Lyapunov / All three actuators 

the beam. This is seen in the response as the actuators, al- 
though stabilizing, drive the beam to excessive amplitudes. 
Small values of k equate to moments of insufficient magnitude 
which decrease the damping ratio. Optimizing the system re- 
quires a tradeoff between the desired settling time and a max- 
imum allowable amplitude. A trial and error procedure is 
required to select the best value and customize the control to 
the application. 

2. Positive Position Feedback (PPF) 

This control algorithm is defined strictly in the 
modal form. For a scalar case, the equations of motion are 
represented by: 
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£ + 2(u>£ + to 2 £ = <?<o 2 ti 



(39a) 



ii + 2C c « c f| + (o 2 c r) = ( 39b > 

where £ is the structure modal coordinate, rj is the compensa- 
tor modal coordinate, and are the structure and compensa- 
tor damping ratios, w and w c are the structure and compensator 
natural frequencies and G is a gain factor. As seen in Equa- 
tion (39) , PPF means positively feeding the structure coordi- 
nate back to the compensator and positively feeding back the 
compensator coordinate times a gain to the structure. 

The system response characteristics of this algorithm 
are shown in Figure 9 . Assuming that the structure motion at 
steady state is of the form: 

£ (fc) = «e iu£ (40) 

the compensator will respond as: 

t| ( t) = pe i(ut ' 40 (41) 

where the phase angle, <£, is defined by: 



4 > 



tan -1 




(42) 
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Figure 9 Positive Position Feedback 

When the structure vibrates at a frequency much lower than the 
compensator natural frequency, the phase angle approaches 
zero. Substituting Equation (41) with <p = 0 back into Equation 
(39a) gives: 

l + 2C«£ + ( (i) 2 - <?Pg> 2 )£ = o (43) 

resulting in a decrease in the stiffness term. When the struc- 
ture and compensator have the same natural frequency, the 
phase angle is 7r/2 . In this case, the structural equation is 
modified to: 

£ + (2Co> + <300)) £ + g> 2 £ = 0 (43) 
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showing an increase in damping. When the structure is vibra- 
ting a frequency much greater than the compensator natural 
frequency, the phase angle approaches tt . Making the same sub- 
stitution yields: 

£ + 2£<i)£ + (to 2 + G0 g> 2 ) £ = 0 (45) 

causing an increase in the stiffness term. These equations 
clearly show that, in order to maximize damping, co c must be as 
closely matched to w as possible. 

This analysis will go beyond the simple scalar case 
and investigate a multi-modal vibration suppression system. 
Specifically, the first three modes will be controlled using 
two different schemes. First, three collocated sensor/actuator 
pairs connected to three compensators tuned to the first three 
modal frequencies will be analyzed. The second system will 
utilize only one collocated sensor/actuator pair connected to 
the same compensators. Figures 10 and 11 illustrate these two 
configurations with simple block diagrams. 

In a multi-modal case, the system equations take the 

form: 

l + Di + Kl = C t GKx\ (4 6a) 

Ti + D c x\ + K c t\ = K c Ci (4 6b) 

For a system with n structural modes and m compensators, £ is 
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Three Actuators and Three Compensators Controlling Three Modes 



SENSORS INPUT OUTPUT 




ACTUATORS 



Figure 10 Three Actuator / Three Compensator Control Scheme 



Single Actuator and Three Compensators Controlling Three Modes 



INPUT OUTPUT 



SENSOR 




ACTUATOR 



Figure 11 Single Actuator/Three Compensator Control Scheme 
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a [n x 1] vector and t| is a [m x 1] vector. D and D c are 
[n x n] and [m x m] diagonal matrices containing the damping 
terms and 2£ ci u> ci . K and K c are [n x n] and [m x m] diag- 

onal matrices containing the squares of the natural frequen- 
cies a> 2 ,- and to 2 ^. G is a [m x m] diagonal matrix of compensator 
gains. C is a fully populated [m x n] participation matrix 
which determines the influence of each sensor/actuator pair on 
each compensator and vice versa. 

The formulation of these matrices is straightforward 
except for the participation matrix. Equating the right hand 
sides of Equations (33) and (46a) : 

& T FU = C T GKx\ (45) 

and assuming the control vector u is equal to the vector GK t] , 
C is equal to (3> T F) T . u is a [p x 1] vector, where p is the 
number of actuators, which forces the dimensions of <i> T F to be 
[n x p] . A modification must be made in the case of the single 
actuator (u = 1) . In order to keep the G and K matrices in the 
same form and ensure dimensional compatibility, a [1 x m] 
vector A whose elements are unity must be included such that: 

<& t Fu = C T AGKr\ (4 8) 
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Now the entire system can be put into the matrix form: 



i 




D O' 


K 




K 


-C t (A) gk 


Y 


A 


+ 


0 D c 


A 


+ 


-K C (A T ) C 




A 



As was done in Chapter II.C.l with Equation (38), a 
state- space transformation is made to perform the simulations. 
The same initial conditions are used to excite the tri -modal 
vibration previously shown in Figure 4. The compensator natur- 
al frequencies were chosen to exactly match the structural 
frequencies in order to maximize the damping effect. In this 
case, the com-pensator damping ratio, f c , should be as small 
as possible. As shown in Figure 12, this creates a resonance 
peak near the compensator frequency allowing the largest gain 
possible to act on the affected mode in the active damping re- 
gion. This will also diminish the effects of active flexibil- 
ity on lower modes. For these simulations, it was assumed that 
each compensator could be designed to have a damping ratio 
equal to .1. In actuality, it will be very unlikely that the 
compensators can be designed as precisely as they have been 
here and a decrease in system performance is expected. To 
prevent a significant decrease, the compensator damping ratio 
should be increased so that the active damping region is wid- 
ened as much as possible. 

Figure 13 shows the response of the beam under the in- 
fluence of three actuators at the base, in the middle and at 
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Figure 12 PPF Compensator Transfer Function Magnitude 

the tip. The effectiveness of this system is seen in how 
quickly all three modes are suppressed. A single actuator at 
the base also provides a sufficient degree of control as seen 
in Figure 14. The effect of actuator placement is evident in 
a comparison with Figures 15 and 16 which correspond to a 
single actuator in the middle and at the tip, respectively, 
feeding the same compensators with the same gains. However, 
the performance of these last two configurations can be im- 
proved by increasing the gains. Due to the fact that a larger 
moment can be placed at the tip than at the base without dri- 
ving the beam unstable, the allowable gains are two orders of 
magnitude greater for Case 3 than Case 1. The individual modal 
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responses can also be tailored by varying the gains such that 
the damping of one mode is maximized in relation to the 
others . 



IIP MOTION - FEM PPF U=3 Gl=.2038 G2=G3-.0i 




TIME, sec 

Figure 13 PPF Three Actuators 

In the scalar case, system stability is clearly de- 
fined. A simple manipulation of the equations of motion and 
application of the Routh-Hurwitz criterion results in the 
condition: 

0 < G < 1 (50) 



This method is not applicable to the much more complex multi- 
modal case. Fanson and Caughey have developed the following 
stability condition: 
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TIP POTION - FEM PPF U=1 e=l Gl=.173 GZ=G3*.BB01 




TIME, sec 

Figure 14 PPF / Actuator at base 



TIP NOTION - FEN PPF U=1 e=5 Gl=.173 G2=G3=.0001 




TINE * sec 

Figure 15 PPF / Actuator in middle 
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TIP rOTlON - FEfl PPT U = 1 e=8 Gl=.173 G2^G3=.0B0i 




TIME, sec 

Figure 16 PPF / Actuator at tip 

K - C t GKC > 0 (51) 

for employing PPF in a multi-modal system [Ref. 8] . This con- 
dition is limited in that this single criterion is a function 
of m variables which are the gains of the individual 
compensators. The procedure used for these simulations was to 
fix two gains while varying the third to determine the range 
where the eigenvalues of the assembled matrix were all 
positive . 

3. Strain Rate Feedback (SRF) 

The final algorithm investigated is also defined in 
modal form as: 
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$ + 2C<d£ + (0 2 £ = -Ga) 2 T] 



(52a) 



n + 2C c w c n + = o> c £ 



(52b) 



where the variables are the same as those previously defined 
for PPF. In this case, the structural velocity coordinate is 
fed back to the compensator and the compensator position 
coordinate is fed back times a negative gain to the structure. 
The gain must be negative, as will be shown in the following 
analysis, so that the system performance will be optimized 
over the most useful frequency range. Figure 17 graphically 
summarizes the performance characteristics of SRF. 

If we again assume that the structural modal coordi- 
nate is of the form: 



£(fc) = oe iut 



(53) 



the compensator response will be: 



i\ (t) = pe i(ot - ♦> 



(54) 



where <f> is now defined by: 




(55) 
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Figure 17 Strain Rate Feedback 

When u << u c , the phase angle, 0, approaches -7r/2. Substitu- 
ting Equation (54) with 0 = -7r/2 back into Equation (52a) 

gives : 

£ + (2C<*> + C?Po) £ + o> 2 £ = 0 (56) 

resulting in an increase in damping. When w = w c , 0=0 and a 
similar substitution yields: 

l + 2 £g>£ + (« 2 + <?p&> 2 )£ = 0 (57) 

increasing the stiffness of the structure. Finally, when u >> 
« c , 0 = tt/ 2 and the structure response is: 
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£ + (2£o> - <?po)) £ + o 2 £ = 0 



(57) 



showing a decrease in damping. The obvious conclusion is to 
design the compensator such that the structural frequencies of 
the modes to be damped fall well below the compensator natural 
frequency. 

The SRF simulations were done using the same format as 
the PPF simulations. The first three structural modes were 
controlled using three compensators and the same four sensor/ 
actuator configurations. Additionally, the possibility of con- 
trolling the first three modes using only one compensator and 
one actuator is investigated. The participation matrix formu- 
lation is identical to the PPF case. The state-space transfor- 
mation is applied to the system defined by: 



r 


+ 


D 0 


T 


+ 


K C t (A)GK 




= 0 


(59) 


A 




-0> c (A r )C D c 


A 




0 K c 









Selecting a precise compensator frequency is not as 
clearly defined as it is with PPF. Additionally, the stability 
condition is not as clearly as defined as in the PPF case due 
to the fact that the closed-loop damping and stiffness matri- 
ces can not be symmetrized. However, the wider active damping 
region gives the designer a lot more flexibility in the cir- 
cuit design. As long as the compensator frequency is greater 
than the structural frequency, a certain amount of active 
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damping will be provided. Although the stabilizing feedback 
gain matrix is heavily dependent upon the structure and com- 
pensator frequencies, some numerical analyses predict that 
higher compensator frequencies produce larger damping ratios. 
A drawback to frequencies in this range is that, as shown in 
Figure 18, the magnitude of the compensator transfer function 
approaches zero limiting the amount of control. Further study 
is required on this point. Again, the compensators were assum- 
ed to have damping ratios equal to .1. To highlight the ef- 
fects of varying compensator frequencies, simulations were run 
for the three actuator/three compensator configuration with 
two different sets of a> c at the two extremes. First, the com- 
pensators were assigned frequencies ten per cent higher than 
each of the structural modes. This response is seen in Figure 
19. Although the higher modes are controlled, the first mode 
is not significantly damped. The system is operating too far 
in the active stiffness region. The dominance of the active 
stiffness can be seen in the higher frequency of the first 
mode. The second set of compensator frequencies were 50 times 
higher than each of the structural modes. This has the effect 
of placing all three structural modes in the active damping 
region of all three compensators. This response is shown in 
Figure 20. In this case, the transfer function magnitudes are 
so small that gains four orders of magnitude larger than the 
previous case are needed. With these large gains, the higher 
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COMPENSATOR TRANSFER FUNCTION MAGNITUDE - STRAIN RATE FEEDBACK 




Figure 18 SRF Compensator Transfer Function Magnitude 



TIP MOTION - TEH SRF U=3 Gl=,0209 G2=G3=.000S 




I IHE ■ sec 

Figure 19 SRF w c = 1.1 x cj 
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TIP MOTION - FErt SRF U=3 Gl=1.6 G2=G3=1 




TlfiEi sec 

Figure 2 0 SRF co c - = 5 0 x co 

modes are almost immediately suppressed. The first mode, hav- 
ing the lowest relative gain, is not controlled as efficient- 
ly. A significant improvement in performance occurs as co c 
moves farther away from to. The same variation in the system 
response when the single actuator moves from the base to the 
tip is seen as in the PPF simulations. 

The control system can be simplified by only using one 
compensator and one actuator at the base as shown in Figure 
21. The motivation for this approach comes from the results of 
the three compensator case where all three structural frequen- 
cies were below the smallest compensator frequency. Theoreti- 
cally, a compensator will provide active damping to all modes 
with a lesser frequency. For this simulation, the compensator 



40 



frequency was set at ten times the third modal frequency. The 
system matrices must be modified because the compensator equa- 
tion is now in a scalar form. As shown in Figure 22, this 
configuration effectively controls all three modes with maxi- 
mum effect on the higher modes. However, this method will only 
work for a few densely spaced modes. If the frequency range is 
expanded too far, the lower modes will not receive a strong 
enough control input. 

This approach was also taken using an alternative 
method. The sign of the gain applied to the compensator ouput 
is changed such that the structural modal equation is now: 

i; + 2 Cg>£ + (i) 2 l = G(o 2 t\ (50) 

This has the effect of swapping the regions of active positive 
damping and active negative damping and making the region cen- 
tered around « c active flexibility. Now, the objective is to 
set the compensator frequency below the first mode placing all 
the modes in the active positive damping region. In this simu- 
lation, the compensator frequency was set at 80% of the first 
modal frequency. A first attempt was made using the same com- 
pensator damping ratio (f c = .1) . This was an unstable system 
due to the resonance peak of the compensator transfer function 
occuring in the active flexibility region. To improve the re- 
sponse, the damping ratio of the compensator had to be in- 
creased such that there was no resonance and the flexibility 
effect was diminished with respect to the damping effect. A 
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SINGLE ACTUATOR AND SINGLE COMPENSATOR CONTROLLING THREE MODES 



SENSOR 



1 

















INPUT 



OUTPUT 



Compensator — , 
SRF 



ACTUATOR 



Figure 21 Single Actuator/Single Compensator Control Scheme 

value for of .95 produced the response in Figure 23. As ex- 
pected, the higher modes are not damped out as quickly as in 
the previous case. To prevent the first mode from being driven 
unstable, the gain must be held relatively low which reduces 
the impact at higher frequencies. In practice, the effective 
bandwidth of the actual compensator must be considered. The 
compensator will damp only those modes densely packed within 
the range of control . 

Deriving a stability condition for SRF similar to 
Equation (51) for PPF proved difficult. A transformation to 
symmetrize the modal equations could not be found. In their 
unsymmetric form, the equations could not be manipulated to 
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Tl? MOTION - FEM SRF U=1 e=l G=.B3 




Figure 22 SRF Single Actuator / Negative Gain 



TIP MOTION - FEM SRF U=1 e=l G=.0B03 




TIME, sec 



Figure 23 SRF Single Compensator / Positive Gain 
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form a simple stability condition. The alternate method used 
was to ensure all the eigenvalues of the system matrix of the 
state- space equation were negative. Again, two gains were held 
fixed while varying the third to determine the stable region. 

4. Comparison of PPF and SRF 

A strict side-by-side comparison of these two methods 
is not feasible due to the fact that the compensators are not 
performing the same functions in both methods. Also, since the 
stable gain regions were so varied, a uniform set of operating 
parameters that provided sufficient control in both systems 
could not be defined. 

Positive Position Feedback is preferable for applica- 
tions where a single mode needs to be quickly damped. Control 
can be maximized by taking advantage of the compensator reso- 
nance occurring in the active damping region. PPF also has the 
advantage of using the sensors to measure the structural dis- 
placement rather than the more difficult measurement of struc- 
tural velocity as is needed with SRF. For multi-modal control, 
the frequency requirements are too restrictive. One compensa- 
tor per mode is required in this case. Unless the structural 
characteristics are uniform and precisely known and unless the 
compensator is closely matched with the modal frequency, PPF 
will be operating at less than optimum. 

Strain Rate Feedback has the advantage of controlling 
modes over a wider range of frequencies. As was demonstrated, 
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a single compensator feeding a single actuator can control 
three modes. This is important in terms of reducing circuit 
complexity and manufacture and implementation costs. However, 
the active damping occurs in a region where the compensator 
gain is not at a maximum. This has the effect of limiting the 
amount of control on modes far from the compensator frequency. 

The robustness of these two control schemes is demon- 
strated in Figures 24 and 25. The scenario simulated here is 
that the actual structural frequencies experienced on-orbit 
are 10% less than the predicted values around which the com- 
pensators were designed. All other parameters are kept con- 
stant. This new response of the "modified structure" is shown 
with the dashed line and is compared to the response of the 
"original structure" shown with the solid line. 

PPF produces almost an identical response for both 
systems with only miniscule differences. Initial expectations 
were that this frequency variation would produce a much more 
significant difference. Further study revealed that the ro- 
bustness of the system is very dependent on the compensator 
damping ratio. Small values produce a phase angle curve that 
approaces a square wave which were seen to emphasize the ad- 
verse effects of a frequency change. Likewise, large values 
stretch out the phase angle curve allowing a less significant 
change in the system response. For this particular model, the 
chosen value is large enough to prevent PPF's frequency depen- 
dency from degrading the system response. The key to providing 
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TIP MOTION - FEM PPF U=3 Gl=.2038 G2=G3=.01 




Figure 24 PPF Structural Frequencies Varied 10% 



TIP MOTION - FEM SRF U=3 Gl=1.6 G2=G3=1 




TIME, sec 

Figure 25 SRF Structural Frequencies Varied 10% 
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sufficient robustness is to have a moderate compensator dam- 
ping ratio. Small values will limit the bandwidth of the com- 
pensator and large values, while still being robust, will 
decrease the magnitude of the damping provided. 

A variation of only the first mode is seen when using 
SRF. The "modified" first medal frequency has moved closer to 
the edge of the bandwidth of the compensators such that the 
amount of damping able to be provided has decreased slightly. 
The structure is also more flexible as seen in the lower fre- 
quency. SRF is much less dependent on the compensator damping 
ratio in terms of robustness. The system is operating in a 
region where the slope of the phase angle curve is not signif- 
icantly altered by variations of the damping ratio. Simula- 
tions showed that the system was equally robust with damping 
ratios different by two orders of magnitude. 
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III. EXPERIMENTAL ANALYSIS 



The characteristics of Positive Position Feedback will be 
studied using the Naval Postgraduate School's Flexible Space- 
craft Simulator (FSS) . 

A. PHYSICAL SET-UP 

The mechanical lay-out of the experiment is shown in Fig- 
ure 26. Reference 4 details the circuit design and the choice 
and placement of the piezoceramic sensors and actuators which 
comprise the active damping control system. A detailed view of 
the sensor and actuator locations is shown in Figure 27. The 
ceramics are mounted on both sides of the flexible arm with 
opposite polarities such that the sensor voltages and actuator 
strains are consistent. The damping circuit consists of three 
sections: a charge amplifier, a low-pass filter and a power 
amplifier. The components of the filter were chosen such that 
the compensator values are w c = .15 Hz and = . 1 . 

For this analysis, the main body of the simulator is held 
fixed and the flexible arm is the only component affected by 
external forces. The flexible arm is supported by air pads 
which allow two-dimensional motion in an approximately fric- 
tionless environment. 
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B . PROCEDURE 



The objective of the experiment was to investigate the 
effectiveness of the PPF algorithm over a range of structural 
frequencies and compare the results with theory. 

To achieve a variety of frequencies, the configuration of 
the flexible arm was changed by attaching and detaching the 
concentrated masses. The natural frequency, to, for a simple 
spring-mass system is defined by the relationship: 



0) 




(61) 



where m is the mass and k is the stiffness. The masses are 
attached to the arm such that the stiffness remains constant. 
Therefore, adding mass will lower the frequency from the ini- 
tial value and vice versa. The selection of where to make mass 
adjustments was a trial and error procedure until frequencies 
of sufficient range and separation were found. 

Data was obtained in the following manner. The signal from 
the sensors was sent via a Gould 5900 Signal Conditioner to a 
Gould RS3800 Stripchart Recorder. This allowed obtaining a 
time history of sufficient length needed to calculate damping 
ratios. Phase angle measurements were made using a Hewlett- 
Packard 54601A Oscilloscope which displayed both the signal 
from the sensors and the signal to the actuators. Frequency 
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measurements were made using both sets of equipment. The data 
was averaged over several run for each configuration. 

With only a single compensator in the feedback circuit, 
the system will only be effective on the mode for which it is 
designed. Therefore, in this experiment, only the first mode 
of the flexible arm was excited so that the cleanest and 
smoothest signal could be recorded. 

The damping ratio, f, was calculated using the log decre- 
ment method: 



C = 



1 

2 iza 




(62) 



where A* is the initial amplitude, A* is the final amplitude 
and n is the number of cycles between those two amplitudes. 
The natural frequency is computed using the damping ratio and 
the damped period, r d : 






27t 



(63) 



However, the damping ratios are small enough with respect to 
the frequency measurement accuracy that we can assume the nat- 
ural frequency to be the same as the damped frequency. The 
phase angle was calculated from the relationship: 

<|> = 6> Afc (3 60° /cycle) (64) 
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where w is the structural frequency in Hz and At is the time 
difference between the sensor and actuator signals as seen on 
the oscilloscope. The theoretical value for the phase angle at 
each structural frequency is calculated using Equation 42. 

C . RESULTS 

As shown in the plot in Figure 28, data was obtained in 
all three regions of active flexibility, damping and stiff- 
ness. (Refer to Figure 9). Although the experimental values 



PHASE ANGLE VS STRUCTURAL FREQUENCY 




- Theoretical 
X Experimental 

Figure 28 FSS Flexible Arm Phase Angles 

for phase angle vary from theory (see Table I) as much as 
seven degrees in the worst case and considering the +/- 2 
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degree measurement accuracy, the data clearly follows the 
trend we expect and we can assume the circuit is functioning 
properly. 

Table I 



FLEXIBLE SPACECRAFT SIMULATOR - PPF 


CASE 


PHASE ANGLE (Degrees) 


Theoretical 


Experimental 


1 


19 


15 


2 


29 


29 


3 


32 


36 


4 


55 


59 


5 


90 


84 


6 


111 


107 


7 


139 


146 


8 


150 


156 ] 



The parameters describing the dynamic response of the 
flexible arm are listed in Table II. Examining the change in 
frequency from the freely vibrating (open-loop) to the active- 
ly damped (closed-loop) cases, the expected effect on the 
stiffness term is evident. At the lower frequencies, a de- 
crease is observed indicating the arm has become more flexi- 
ble. At and near the compensator frequency, there is little or 
no change. The higher frequencies exhibit an increase meaning 
the arm has become stiffer. 
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Table II 



FLEXIBLE SPACECRAFT SIMULATOR - PPF 


CASE 


FREQUENCY (Hz) 


DAMPING RATIO 


Free 


Damped 


Free 


Damped 


% A 


1 


. 116 


.113 


.0095 


.0070 


-27 


2 


. 128 


.125 


. 0104 


. 0114 


10 


3 


.133 


.128 


. 0077 


.0243 


215 


4 


.141 


. 140 


. 0059 


. 0325 


450 


5 


. 150 


.150 


.0059 


.0221 


345 


6 


. 154 


.156 


.0040 


.0158 


299 


7 


. 165 


.168 


. 005? 


. 0145 


161 


8 


. 177 


. 180 


. 0035 


. 0093 


165 | 



A very pronounced variation is seen in the change in the 
damping ratio when the feedback system is activated. At the 
lowest observed frequency, the arm has become so flexible that 
the damping worsens by 27 per cent. There is a steady improve- 
ment as the structure moves into the active damping region. 
The maximum effect is observed with Case 4 showing a 450 per 
cent improvement . The magnitude of the control provided is 
best appreciated visually as seen in Figure 29. Even at the 
higher frequencies, while not in the active damping region, an 
improvement in the damping ratio is produced by the increased 
stiffness of the arm. 
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One discrepancy evident in this data is that the greatest 
improvement in the damping provided by the compensator does 
not occur at a structural frequency equal to the compensator 
frequency. The cause of this variation has been difficult to 
pin-point. In order for the transfer function resonance peak 
to occur at .141 Hz, where the maximum increase was observed, 
the compensator damping ratio would have to increase to .34 or 
the compensator frequency would have to decrease to .142 Hz. 
However, this is not consistent with the phase angle data pre- 
viously mentioned. This rules out any small variations in the 
feedback circuit as possible sources of error. There are some 
physical factors which may have caused inconsistencies between 
the eight configurations. Other experiments were being conduc- 
ted while this one was in progress which did not allow a 
uniform orientation of the flexible arm for all of the config- 
urations. Varying frictional or gravitational effects may have 
adversely affected the consistency of the results. Further de- 
tailed experimental work is needed to precisely determine 
where the maximum effect actually occurs. 

This data highlights the advantages and disadvantages of 
the Positive Position Feedback algorithm. When the compensator 
is precisely tuned to the structural frequency, the system is 
operating near the maximum gain of the compensator transfer 
function. This enables the active damping control to have its 
greatest effect. However, any deviation from this optimum re- 
gion reduces the impact of the control system. As seen in this 
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experiment, only a 25 per cent reduction in the structural 
frequency results in the compensator degrading the dynamic re- 
sponse of the system. If the frequencies can not be matched 
precisely, it is better to err into the active stiffness re- 
gion where the structural characteristics are not adversely 
affected. The efficiency of the control law with a bounded 
gain space, 0<G<1, is heavily dependent upon the compensator 
parameters <o c and f c . It is recommended that moderately higher 
compensator damping ratios be used so that a relatively large 
amount of damping is provided while still maintaining suffi- 
cient robustness. 



57 



IV. CONCLUSIONS 



Piezoceramic sensors and actuators have been shown to pro- 
vide an effective means of controlling structural vibrations. 
An increase of 450 per cent in the damping ratio of the flexi- 
ble arm of the Flexible Spacecraft Simulator was demonstrated 
using the Positive Position Feedback algorithm. This system 
provides the greatest degree of control when the compensator 
is precisely tuned to the modal frequency. Its main limitation 
is that only a single mode can be controlled per compensator. 

Positive Position Feedback and Strain Rate Feedback have 
significantly different control characteristics which were 
analytically demonstrated for a cantilevered beam. SRF is more 
advantageous than PPF in the multi-modal environment. With a 
wider active damping region, it has the ability to stabilize 
more than one mode given a sufficient bandwidth. 

Recommendations for further study include a number of 
topics. First, the dynamic response of the FSS flexible arm 
could be improved by incorporating additional compensators 
into the current configuration. The first mode is very rarely 
the only mode excited during slew maneuvers and the higher 
modes should be considered. This problem may be more easily 
solved by building a Strain Rate Feedback circuit and integra- 
ting it with the sensors and actuators currently mounted on 
the arm. Also, the effect of multiple sensors and actuators on 
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the arm can be investigated with an optimization study of sen- 
sor/actuator placement. There are a variety of control laws 
which can be used in this application. As an example, a combi- 
nation of PPF and SRF could be developed which may incorporate 
the strengths of each and cancel out the weaknesses. 

The Flexible Spacecraft Simulator is an excellent tool for 
studying, not only active damping control techniques, but 
other spacecraft attitude control and dynamics problems. Cur- 
rent research endeavors include designing an improved thruster 
for attitude control and incorporating a robotic arm onto the 
main body. 
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APPENDIX 



% Simulation of Beam Vibration using Finite Element Method 
% Transformation into modal form using PPF with three control 
laws 

% utilizing three actuators on elements 1,5,8. 

% 

% Define physical constants: 

% 



n=8 ; 
h= . 1 ; 
d= . 5 ; 
EI= . 9 ; 
inertia 



% number of elements 
% element length 
% element density 

% modulus of elasticity times moment of 




% Element mass matrix: 



% 

Mul= (d*h/420) * [156 22*h; 22*h 4*h A 2] ; 
Mur= (d*h/420 ) * [54 -I3*h; I3*h -3*h A 2] ; 
Mll=Mur' ; 

Mlr= (d*h/420) * [156 -22*h; -22*h 4*h A 2] ; 

% 

% Element stiffness matrix: 



% 

Kul= (EI/h"3) * [12 6*h; 6*h 4*h A 2]; 

Kur= (EI/h A 3) * [-12 6*h; -6*h 2*h A 2] ; 

Kll=Kur ' ; 

Klr= (El/h A 3) * [12 -6*h; -6*h 4*h A 2] ; 

o. 

*o 

% Combine matrices for n element system. Disregard first two 
rows 

% and columns for clamped end. 

% 

for i=l:n-l, 

M(2*i-l:2*i,2*i-l:2*i) =Mul+Mlr; 

K(2*i-l:2*i,2*i-l:2*i) =Kul+Klr; 

% 



M(2*i-l:2*i,2*i+l:2*i+2) =Mur; 
K(2*i-l:2*i,2*i+l:2*i+2)-Kur; 

% 

M(2*i+l:2*i+2, 2*i-l:2*i) -Mil; 
K (2*i + l:2*i+2,2*i-l:2*i) =K11; 

end 



M(2*n-l:2*n,2*n-l:2*n) =Mlr; 
K(2*n-l:2*n,2*n-l:2*n) =Klr; 

% 

Mi=inv(M); % Inverse of mass matrix. 
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% Obtain eigenvalues and eigenvectors: 

% 

[phi, lambda] =eig (K,M) ; 

% 

%• Extract natural freqencies: 

% 

w=diag (lambda) ; 

% 

% Sort from low to high frequency: 

% 

[wl , num] =sort (w) ; 
w=w (num) ; 
for j =1 : 2 *n 

lambda ( j , j ) =w ( j ) ; 

end 

% 

% Sort eigenvectors: 

% 

phi=phi ( : , num) ; 

% 

% Calculate generalized masses: 

% 

Mg=phi ' *M*phi ; 

% 

% Calculate weighted modal matrix: 

% 

for i=l:2*n 

phiw ( : , i) =phi ( : , i) /sqrt (Mg(i, i) ) ; 

end 

Mw=phiw' *M*phiw; 

Kw=phiw' *K*phiw; 

% 

% Actuator location matrix (actuators on elements 1,5,8) 



F=zeros (16,3) ; 

F(2,1)=1;F(8,2)=-1;F(10,2) =1;F (14,3) =-l;F (16,3) =1; 

Ft=phiw' *F; 

% 

% Assign structural modal damping values: 

% 

Z (1) = . 0001; z (2) =. 0004; z (3) =. 0009; z (4) = . 0015; z (5) = . 0037; Z (6) = 
.0089; 

z(7)=.01;z(8)=.025;z(9)=.042;z(10)=.067;z(ll)=.098;z(12)=.ll; 
z (13) = .35; z (14) = .86;z (15) =1.21; z (16) =1.77; 

% 

% Assemble structural damping matrix: 

% 

D=zeros (2*n) ; 
for i=l:2*n 

D(i,i)=2*z(i) *sqrt (Kw(i, i) ) ; 
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end 

% 

% Truncate to first three modes: 

% 

Kwl=Kw(l:3, 1:3) ; 

Ftl-Ft (1:3, 1:3); 

D1=D ( 1 : 3 , 1 : 3 ) ; 

% 

% Assign compensator values: 

% 

wcl=sqrt (Kwl (1,1) ) ;wc2=sqrt (Kwl (2,2) ) ;wc3=sqrt (Kwl (3,3)); 
zcl= . 1 ; zc2= . 1 ; zc3= . 1; 

Gl= . 2038 ;G2= . 01 ; G3= . 01 ; 

Kc= [wcl^2 0 0; 0 wc2 > '2 0; 0 0 wc3^2] ; 

Dc= [2*zcl*wcl 0 0; 0 2*zc2*wc2 0; 0 0 2*zc3*wc3] ; 

G=[G1 0 0; 0 G2 0; 0 0 G3] ; 

% 

% Combine structure and compensator equations into matrix 
form: 

% 

Ds= [D1 zeros (3); zeros (3) Dc] ; 

Ks= (Kwl - Ft l*G*Kwl ; -Kc*Ftl' Kc] ; 

o. 

*© 

% Put closed- loop system in state space form: 

% 

A= [zeros (6) eye (6); -Ks -Ds] ; 

B=zeros (12,1) ; 

C=eye (12) ; 

D=zeros (12,1) ; 

T= (0: .01:9 .99) ; 

U=zeros (T) ; 

X0=zeros (2*n, 1) ; 

% 

% Initial displacements to excite 3 modes 
% 

X0 ( 1 , 1) = . 008 ;X0 (2 , 1) = . 05 ;X0 (3 , 1) = . 01 ;X0 (4 , 1) =0 ;X0 (5 , 1) = . 004 ; 
X0 (6, 1) =- . 086 ;X0 (7, 1) =- .007;X0 (8, 1) =- .086;X0 (9, 1) =- .016;X0 (1 
0 , 1 ) = 0 ; 

X0 (11, 1) =- .015;X0 (12, 1) =.01;X0 (13, 1) =- .001;X0 (14, 1) =.173; 

X0 (15,1)-.02;X0 (16, 1)-.173; 

Z0=inv (phiw) *X0 ; % Initial displacement in modal 

coordinates 

Z01=Z0 (1:3) ; 

M0= [Z01' 00000000 0]'; 

[Y, X] = 1 s im (A,B,C,D,U,T, M0 ) ; 

% 

% Put open- loop system in state -space form: 

% 

Aol= [zeros ( 3 ) eye (3); -Kwl -Dl] ; 

Bol=zeros (6,1) ; 

Col=eye ( 6 ) ; 
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Dol=zeros (6,1) ; 

MOol= [Z01 ' OOO]; 

[Yol , Xol] =lsim (Aol , Bol , Col , Dol , U, T, MOol) ; 

% 

% 

% Convert back to tip physical coordinate: 

% 

phiq=phi w (1:16,1:3) ; 
for i=l:1000 

tm= [Y(i, 1) Y ( i , 2 ) Y(i,3)]'; 

q=phiq*tm; 

tip (i) =q ( 15 ) ; 

tmol= [Yol (i, 1) Yol(i,2) Yol(i,3)]'; 
qol=phiq*tmol ; 
tipol (i) =qol (15) ; 

end 

% 

% Plot data 
% 

plot (T, tip ( : ) ) 
xlabel ('TIME, sec') 
ylabel (' TIP DISPLACEMENT' ) 

title ('TIP MOTION - FEM PPF U=3 Gl=.2038 G2=G3=.01') 
pause 

plot (T, tipol ( : ) ) 

xlabel ('TIME, sec') 

ylabel (' TIP DISPLACEMENT' ) 

title ('TIP MOTION - FEM OPEN LOOP') 

pause 
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